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SUMMARY 


A method which provides sonic-boom-minimizing equivalent area distributions 
for supersonic cruise conditions is described* This work extends previous anal- 
yses to permit relaxation of the extreme nose bluntness which is required by 
conventional low-boom shapes and includes propagation in a real atmosphere* 

The procedure provides area distributions which minimize either shock strength 
or overpressure. The minimizing Whitham F-f unctions for the pressure signatures 
are also included. 


INTRODUCTION 

Sonic-boom levels produced by future generation supersonic transports must 
be considered if their supersonic flight paths, and thus their economic success, 
are not to be severely restricted. It is generally believed that the best 
approach to this problem is to include low-boom constraints early in the design 
process. The procedure described herein provides this constraint in the form 
of an equivalent area distribution which depends on the design conditions of 
aircraft weight (lift) and length, and on cruise Mach number and altitude - 

Sonic-boom-minimization studies, by Jones, Seebass, and George for uniform 
and isothermal atmospheres have been reported in references 1 to 6* These anal- 
yses have been extended to a real atmosphere in reference 7. All these studies 
have indicated that the minimizing area distribution is characterized by a blunt- 
ness in the leading portion of the area distribution, i.e., the aircraft nose 
region. Because extreme nose bluntness produces large drag, a method of relaxing 
the bluntness requirement is needed to offer the opportunity for compromise between 
blunt-nose, low-boom and sharp-nose, low-drag configurations. 

This analysis extends the work of Seebass and George (ref. 6) and Darden 
(ref. 7) to include nose-bluntness relaxation. Because questions still remain 
regarding the type of pressure signature and the level of overpressure which 
are acceptable, either a minimum-shock or a minimum-overpressure signature may 
be produced for given cruise conditions and nose length. 

This analysis has been incorporated into a computer program which is 
described in the appendix of this paper. 


SYMBOLS 

Although values are given in both SI and U.S. Customary Units, the calcu- 
lations for the investigation were made in U.S. Customary Units. 

Notation used in the computer printout is given in parentheses. 


A 


Ae(n 

a 

B 


C 

Cd 

D 

F (F and 

G 

H 


h (Z) 
I 


ray-tube area 
equivalent area 

equivalent area at equivalent length I 
speed of sound 

slope of rise in F-function (eq. (20)) 
height of F-function at balancing point 
drag coefficient 
constant in F-function equation 
FTAU) Whitham F-function 
area of F-function to front balance point 
height of spike in F-function 
airplane altitude 
impulse of pressure signature. 


'^P>0 


P dt 


K (RFK) reflection factor 

I (L) equivalent length of airplane 

M Mach number 

P pressure perturbation 

p ambient pressure 

Ap (Delta P) overpressure 

S slope of balancing line 

t time 

u speed 

V total area, gW/pu^ 

W (WG) airplane weight 

V 

w (W-TILDE) weight parameter, - — 

B(1 - yf)^/2 


2 


x,y axial distances 

Yf (YF) width of forward spike in F-function 

Yj. (T) position of rear area balancing 
2 (r) vertical distance of signal from airplane axis 
ot advance of acoustic rays 

6 = (m2 - 1 ) V2 


6 z (BR) 
Y 

r 

6 

n 

X (LAM) 

y 


shifting factor for axial position 
ratio of specific heats, 1 .4 for air 

= ^ ^ ^ 

2 

Dirac delta function 

fraction of balancing line slope used in defining F-function 
axial position at which F-function becomes negative 

Mach angle, sin ' - 

M 


K dummy integration variable 

P density 

Subscripts: 
a ambient 

f front shock 

g ground 

h altitude of initial waveform 

o minimum-overpressure signature 

r rear shock 

s minimum-shock signature 


y 


axial position 



yf axial position of front area^ balance 

z vertical distance of signal from airplane axis 

First and second derivatives with respect to distance are denoted by single 
and double primes, respectively. 


BACKGROUND 

Sonic-boom-miniraization studies are based primarily on currently accepted 
sonic-boom prediction methods. These methods, outlined in figure 1, are exten- 
sions of work done by Whitham (ref. 8), Walkden (ref. 9), and Hayes (ref. 10). 
Prom the complex aircraft configuration and its lift distribution, an equivalent 
area distribution is defined. This equivalent area distribution then defines 
the Whitham "F-f unction" through the following integral relation; 


F(y) 


1 py 

2tt Jq 


fte" 

(y - a'/' 


<3^ 


(1) 


This F-function represents a distribution of sources which causes the same dis- 
turbances as the aircraft at some distance from the aircraft. 

Because the pressure signal propagates at the local speed of sound, and 
each point of the signal advances according to its amplitude, the signal is dis- 
torted at the ground and, theoretically, can be multivalued. The physically 
unrealistic multiple values of pressure are, however, eliminated by the intro- 
duction of shocks. Shock location, based on the observation that for weak dis- 
turbances the shock bisects the angle between two merging characteristic lines, 
is determined by a balancing of the signature areas within loops on either side 
of the shocks. This procedure is demonstrated by the shaded areas of the dis- 
torted signal in figure 1 . 

For present-day supersonic aircraft, the propagation distance (altitude) 
and the pattern of shock coalescence have been such that, at ground level, only 
two shocks remain in the signature. These two shocks have a linear pressure 
variation between them (as seen in fig. 2), thus the name "far-field N-wave." 

For this wave, the shape of the generating aircraft has an effect on the magni- 
tude but no effect on the shape of the resulting signature. For a sufficiently 
long and slender aircraft, it was found that the ground wave form may possibly 
not have attained its N-shape remaining instead a mid-field wave. Because the 
shape of the aircraft does influence the shape of a mid-field pressure signature, 
aircraft shaping has now become a powerful tool in reducing the sonic boom. 

The F-f unctions for the lower bound of an N-wave (refs. 1 and 2) and for 
the lower bound of the bow shock in a mid— field signature (refs. 3 and 4) led 
to the form of the minimizing F-function assumed by Seebass and George for the 
entire signature. Lung (ref. 11) later proved this to be the minimizing form 
by "bang-bang" control theory. 
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Common to all of these minimizing forms of the F-function is a Dirac delta 
function at y = 0. This infinite impulse corresponds to an infinite gradient 
at the nose of the equivalent area distribution. Aircraft designed to match 
these area distributions generally have extremely blunt nose shapes which result 
in substantial drag. This result, though seemingly paradoxical, can be explained 
by the created shock-attenuation pattern in which the shock strengths, and there- 
fore the drag, are greatest near the aircraft. Because of special shaping and 
area growth, secondary shocks from other aircraft components do not overtake and 
enhance the bow shock during the propagation of the wave front. (See fig. 3.) 

The net result of this process is weaker shocks at larger distances because of 
attenuation, but an overall increase in drag. In contrast, notice that the bow 
shock is weaker on the sharp-nosed aircraft, but coalescence with stronger shocks 
from rearward components causes a much stronger shock at mid- and far-field dis- 
tances. If an aircraft configuration is to incorporate boom-minimization features 
and not suffer aerodynamic penalties which inordinately affect its flight effi- 
ciency, range, payload, productivity, and so forth, then means must be provided 
for including nose-bluntness sensitivity in design trade-off studies. 


ANALYSIS 

The following is a brief description of the Seebass-George minimization 
scheme with modifications to provide propagation through the real atmosphere 
(ref. 7) and control over the bluntness of the area distribution and thereby 
the drag of the configuration. The assumed form for the class of minimizing 
F-functions is shown in figure 4. In mathematical terms it may be expressed 
as: 


F(y) = 2yH/yf 

(0 ^ y ^ yf/2) 

(2a) 

F(y) = C(2y/yf - 1) - H(2y/yf - 2) 

(yf/2 = y < yf) 

(2b) 

F(y) = B(y - y^) + C 

(yf = y < X ) 

(2c) 

F(y) = B(y - yf) - D 

(X $ y < l) 

(2d) 


In these equations H, B, C, D, and X are unknown coefficients which are 
determined by the cruise conditions of the aircraft, by nose length, by the pre- 
scribed ratio of bow to rear shock, and by the signature parameter being mini- 
mized. Types of signatures studied include flat-topped signatures in which over- 
pressure is minimized with B = 0 and signatures in which F(y) is allowed to 
rise between yf and X with a resulting minimum shock followed by a pressure 
rise as illustrated in figure 4. The value of B in this form of F(y) may 
range between 0 and S. Recall that the Whitham function F(y) represents the 
shape characteristics of the pressure signature and is defined in reference 10 
and equation (1) in terms of the equivalent area distribution as: 


F(y) 



(y - 
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Crucial to this minimization technique is the fact that equation (1) is an Abel 
integral equation which may be inverted to give the function Ag in terms of 
the F-function. When this function is evaluated at 2, the result is: 

Ag(2) = 4 r F(y)(2 - y) ^ /2 <3y (3) 

Upon substituting the minimizing form of the F-function into equation (3) and 
integrating, the following equation for the development of cross-sectional area 
is obtained; 


Ap(x) 




3/2 


15 yf 


2/15' 


'3yf \ 

— + j(2C ~ 4H) 


+ 5(2H - C) 


+ 1 (X - yf) 4(x - yf) 


3/2 


2c\f 2 \ 2 

.-c 


4 /H \ 4 2 2 2 

+ — — (3yf + 2x) - - H + — B(3yf + 2x) - - Byf + - C 

15\yf/ 3 15 33 


1 (X - X)-(x - X)3/2(c + D) 
3 


(4) 


where 1 (x - 2) is the Heaviside unit step function. A typical form of the 
resulting area distribution, Ag, is seen in figure 4. 

If the effects of aircraft wake and engine exhaust are neglected, and, if 
the aircraft cross-section area is zero at its base, then the area at 2 is 
entirely due to cruise lift or 


W 32 H 8 / 

A ( 2 ) = 2 5/2 + _ h 

® pu2 1 5 yf 1 5 \ 2 ) 


'3yf\ /l \ 

+ 22 — (2C - 4H) + 5(2H - C) 

. 2 / \yf / 


+ 4(2 - yf) 


3/2 


-V- . 2!, + i c + * 2!) - ' B 


.yf 


3 15\yf/ 


2 2 2 

+ — B(3yf +22) Byf + - C 

15 3 3 


_ - X)3/2jc + D) 
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(5) 
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The first constraint imposed upon F(y) is that the front area balance must 
occur at y = yf, where yf is the first point at which F(y) = C, that is, 


pyf ayf ayf 

\ F(y) dy = G = F(yf) = C (6) 

Jq 2 2 


For the real atmosphere, the advance (ref. 9) of any point of the signal is given 
by 


F(y) j^z ^pah\V2/Aj^ ^/2 „ 

(2B)^/2 Jq Apha/ V^hV B 


where z is the vertical distance of the signal from the aircraft axis, 

is the vertical distance of the initial waveform from the aircraft axis, A and 

Ajj are ray-tube areas determined from 



( 8 ) 


and r = (y T)/2. The initial waveform must be defined away from the aircraft, 
because the F-function only can represent the body shape accurately at several 
body lengths away, and the acoustical theory used in describing the propagation 
of the signal fails near the aircraft. The slope of the balancing line is pro- 
portional to the reciprocal of the advance at any point of the signal; thus. 


(26) V2 



F(y) 

ay 


(9) 


For the rear area balancing to occur between points l and then 




yr 


F(y) 



- yf) - D + F(y^)] (y^ - 2) 


( 10 ) 


where yj^ is the unknown second intersection point of the rear area balancing 
line with P(y) . If a cylindrical wake is assumed, F(y) and its integral for 
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y > I can be expressed in terms of F(y) for y < I according to reference 1 
as follows: 


F(y) 


1 


7T(y - 2)^/2 



(Z - C)V2 
y - C 


F(^) 


(y > Z) (11) 


pyr 2 Jyr - 

\ F(y) dy = - ^ F(Q j 


(y > n (12) 


It is necessary to define F(y) and its integral for y ^ 1 in this way for 
optimization problems, since aircraft geometry and thus the Whitham function 
can be varied arbitrarily only in the range 0 ^ y ^ I. The constraint on the 
ratio of shocks is given by 


Pj. D - B(2 - Yf) + F(yj-) 

To insure that y^ is an intersection point of F(y) and the balancing line, 
then 


F(yr) = S(yr - 2) + B(2 - y^) - d 


( 14 ) 


and the slope of F(y) at y^ must be less than S. 

Solving the system of equations ((5), (6), (10), (13), and (14)) provides 
values for the constants H, C, D, X, and y^. One difficulty of this mini- 
mization approach is that, as yet, there is no precise definition of the rise 
time that will allow the ear to detect only the bow shock and not to be sensi- 
tive to the peak as well. With the values of the coefficients known, the mini- 
mizing F-f unction and the area distributions may be determined by equations (2) 
and (4) . 

To convert F(y) into pressure near the airplane, the following equation 
is used: 



( 15 ) 


and on the ground 
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( 16 ) 


P 


g 



Pgag\V2 

Ph 

Ph^h/ 


Finallyr by using a ground reflection factor K, the pressure perturbations are 
converted into the ground signature by 

Apg = PgK (17) 

The ground-level signature is presented by the signature shown in figure 4, 

By taking the limit as yf 0 and by using L'Hospital’s rule, the expres- 
sion for area development given above (eq. (4)“) reduces to the following delta 
function form as given by Seebass and George (ref. 6) : 


A_(x) = 4 GxV2 + — Bx 5/2 + - 0x^/2 - i (x - X)-(x - X)V2(c + D) (18) 

15 3 3 


SPECIAL CASES 

Special types of pressure signatures occur for large or small values of 

~ Bw 1 

the weight parameter w, where w = ■ ■ — ■ . These special cases 

pu^ B(2 - Yf) 5/2 

include signatures with no bow shock, signatures with no shocks, and signatures 
in which the rear shock is less than the bow shock without restriction (fig. 5). 
Values of w may be checked to indicate when the special calculations are 
necessary. 


No Bow Shock 


For no bow shock, the expression for total area r'educes to the yf = 0 
form (eq. (18)) automatically. Using the facts that G = ac2 and X = Z at 
the minimum length necessary for no shock, the resulting quadratic in C may 
be solved explicitly to give 


C 



— 

15 


n 

+ - 


V 


2 B2 5/2, 


1/2 


where 


V = 



B = Sri 



( 19 ) 


(0 ^ n S 1 ) (20) 
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For no bow shock (C = 0) equation (1 9) reduces to 


V 16 

w = i — 

Bj5/2 15 


( 21 ) 


in the minimum shock signature. 


No Shocks 

To eliminate both shock waves, G and C are zero and there must be no 
discontinuity in F at 2. In fact, both F(y) and F’ (y) must be continuous 
at 2 and thus Ag’(x) and Ag"(x) must be zero there, since Ag is constant 
for X > 2. With no shocks. 


16 _ 8 _ 

Ag(2) = — 32^/2 - - D{2 - X)3/2 (22) 

® 15 3 


By taking the first and second derivatives and solving the two resulting equa- 
tions simultaneously 


X 



(23) 


2n 

D = — 
3a 


(24) 


Substituting these values back into the expression for Ag gives 


V 

< 0.47407 

B2 5/2 


(25) 


No Restriction on Rear Shock 

Values of X < 2 in the minimizing F-f unction result because of restric- 
tions on the rear shock. For short lengths, high weights, etc., or other combina- 
tions which give large values, of w, the rear shock is less than the bow shock 
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without restriction. In these instances X = 2 and the system of equations is 
solved without restricting the ratio of rear to front shock; 


2yh 

F(y) = 

Yf 



F(y) = B(y - yf) + C 


(0 < y S yf/2) 

(26a) 

(Yf/2 < y < yf) 

(26b) 

(yf/2 < y < I ) 

(26c) 


Substituting this expression for F into equation (3) and integrating gives a 
quadratic equation in C which can be solved explicitly. Values for C, h, 
X, 2, and B are now established. For the rear area balance to occur at y^ 



F(y) dy = ~[b( 2 - yf) 


D + F(yr)] (yr - 2) 


(27) 


also, yj. must be an intersection point for F(y) and balance point, or 


F(yr) = s(yj- - 2) + B(2 - yf) - D 


(28) 


Combining these equations 



(29) 


which may be solved iteratively for yj. with 


F(yr) 


TT(yr 


2)V2 



(2 - pV2 

— F(^) d^ 

Yr - ^ 


(30) 
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and 


I" 


F(y) dy 


2 ,/yr - 

= - - F(^) tan“M : T-) <3^ 

IT <4l ' 




(31) 


Substituting back into equation (28) gives the required value for D, 

Values of w producing the special cases are shown in figure 5. Typical 
values giving these extreme w's at M = 2.7, z = 18 288 m (60 000 ft) , 

W = 272 155 kg (600 000 lb) would be I = 33.5 m (110 ft) for no restriction 
on the rear shock and I = 167.64 m (550 ft) for no bow shock. These values 
may vary slightly with yf as seen. 


TYPICAL RESULTS 

Except as indicated otherwise, all results shown are for the following con- 
ditions: Mach number, 2.7; altitude, 18 288 m (60 000 ft); weight, 272 155 kg 

(600 000 lb); equivalent length, 91.44 m (300 ft); reflection factor, 2.0; ratio 
of bow to rear shock pressure, 1; and B = 0.5S. The value of the front balancing 
point yf determines the length of the conical nose region of the equivalent area 
distribution and is the controlling factor for drag when all of the other flight 
conditions remain constant. Referring again to figure 4, this position is indi- 
cated on the F-f unction and the corresponding area distribution. 

The effect that varying yf has on the equivalent area distributions for 
both the minimum-shock and minimum-overpressure signatures is seen in figure 6. 
Each of the distributions was determined from the inversion formula which defines 
the effective area corresponding to the minimizing F-function for a given value 
of yf. Values of the overpressure and impulse which correspond to these distri- 
butions have been inserted for reference. Recall that the impulse is the integral 
of the positive portion of the pressure signature, as indicated in the signatures 
at the top of figure 6. Note that the lowest values of overpressure and impulse 
for these conditions occur when the area distribution has an infinite gradient at 
the nose. As yf increases, the longer conical nose region, generated to reduce 
drag, also causes an increase in overpressure and impulse. Thus, extensive trade- 
off studies between drag and boom levels would be necessary in any aircraft design 
studies. The ratio of the increase in overpressure and impulse for both types of 
signatures as a function of the ratio of nose length to airplane length yf/1 is 
seen more easily in figure 7. With increasing values of yf, note that the area 
under the peak portion of the F-function also increases but not in the same pro- 
portions. Therefore, to achieve the necessary total area, the level of the con- 
stant portion of the F-function also must increase, thereby producing a higher 
level of Ap and I. For illustrative purposes, the nose spike of the F-function 
for yf = 0 has been drawn with a finite width. For this condition, the spike is 
defined mathematically as a pulse of zero width and infinite height which never- 
theless has a finite area. In these figures, Ap is the level of overpressure at 
the bow shock and represents the initial "bang" heard by the ear when a sonic boom 
occurs. 
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The variations of overpressure level and the corresponding drag levels are 
seen in figure 8. To develop the corresponding variation of drag with yf, the 
assumption is made that necessary configuration changes would be confined to 
the fuselage forebody itself. The increment in drag produced by these changes 
was calculated using the near-field program of reference 12 and then applied to 
drag values for a typical arrow-wing SST configuration. As seen in figure 8, 
there is a sharp decline in drag as the blunt nose is converted to a conical 
region. Thus, the boom levels could be reduced significantly without prohibi- 
tive drag penalties by defining the proper ratio yf/2. 

The actual trend of overpressure with length is seen in figure 9. The low- 
est curve represents the yf = 0 case, the higher ones increasing with yf. 

In all instances observe that overpressure levels decrease with length. Though 
these results are for minimum-overpressure signatures, similar trends are found 
to exist for the minimum-shock signatures. 


CONCLUDING REMARKS 

The analysis for a sonic-boom minimization method has been presented. The 
method, which represents an extension of previous work, provides the minimizing 
equivalent area distribution with relaxation of the nose-bluntness requirement 
for given cruise conditions and includes propagation in the standard atmosphere. 
The minimizing Whitham F-functions for the minimum-overpressure or minimum-shock 
signature are also included. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
November 13, 1 978 
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APPENDIX 


COMPUTER PROGRAM DESCRIPTION AND FLOW CHART 

This program calculates the Whitham F-f unction and the corresponding eguiv- 
alent area distribution which produce the minimum-overpressure or minimum-shock 
signature* The solution is based on five governing eguations: total area growth 

(eq. (5)), front area balancing point (eq* (6)), rear area balancing point 
(eq* (TO)), ratio of front to rear shock (eq. (13)), and the area balancing line 
at the rear balancing point (eq. (14)). These five equations reduce to two non- 
linear equations in two unknowns which are solved iteratively using the Newton- 
Raphson and a combination of the secant and bisection methods. 

In calculating the advance of the wave front the standard atmosphere has 
been used except for an area near the aircraft axis where it is necessary to 
assume a uniform atmosphere. 

This program was coded for the Control Data 6600 computer in Fortran IV. 

The field length required is 100 000 octal units. Program SEEB is available 
as LAR 11979 from 


COSMIC 

Suite 112 Barrow Hall 
University of Georgia 
Athens, GA 30602 

The flow chart of program SEEB is on the following page. 


INPUT DESCRIPTION 
Cruise Mach number (1.25 to 3.0) 

Cruise altitude (20 000 to 100 000 ft (6096 to 30 480 m) ) 

Airplane equivalent length (200 to 700 ft (60.96 to 213.36 m) ) 

Gross weight of airplane (200 000 to 1 000 000 lb (90 720 to 
453 600 kg)) 

Balancing point in F-function to determine the front shock 
(0 ^ yf/1 ^ 0.2) (Default =0.0 delta function) 

= 0 Iterations on Lambda and T not printed. 

= 1 Iterations printed. 

(Default = 0) 

RKF Reflection factor. 

(Default = 2.0.) 


M 

Z 

L 

WG 

YF 

I PRINT 
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APPENDIX 


BLEN 

PRPF 

WIWG 

ISIG 

IS 

PERCEN 

ICH 

ILEN 

lYF 


Factor determining where iteration on T begins. T = BLEN*L 
(Default = 3.0) 

Ratio of rear to front shock 
(Default = 1 .0) 

Percentage of gross weight remaining at start of cruise. 

W = WIWG-WG (Default =1.0) 

= 1 Area, F-f unction and pressure signature are printed. 

= 0 Only the pressure signature is printed. 

(Default = 0) 

= 0 Minimum-shock and minimum-overpressure signatures. 

= 1 Minimum-overpressure signature only. 

= 2 Minimum-shock signature only. 

(Default = 0) 

Factor determining the difference in the slope of the balancing line 
and the rise portion of the minimum-shock F-function. Slope of rise 
= Slope of balancing line minus PERCEN times slope of balancing line. 
(Default = 0.5) 

Indicates whether the altitude or Mach number changed after the initial 
run and the advance must be recalculated. 

= 0 No change. 

= 1 Mach number or altitude did change. 

(Default = 0) 

= 0 Find Ap for given input conditions - 

= 1 Iterate to find lengths necessary to give Ap = 0.5, 0.6, 0.7, 

0.8, 0.9, and 1.0 psf. A full set of flight conditions includ- 
ing length must be input. 

(Default = 0) 

= 0 Delta function used in F-function. 

= 1 Nonzero value of y^ used. 

If the value of yf is less than 1 , the program automatically uses 
the delta function in the F-function. 

(Default = 1 ) 
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OUTPUT DESCRIPTION 

The output consists of a listing of input conditions, the solution coeffi- 
cients, solution tables, and summary values. Codes exist so that only the pres- 
sure signature is printed, the F-function, area, and pressures are printed, or 
extensive printout occurs during iterations. 

The axial position of the front balancing point is yf. 

The advance factor on the ground is ot. 

Altitude is given in feet and kilometers. 

Gross weight is given in pounds. 

Cruise weight is given in pounds and kilograms. 

WIWG = Ratio of beginning cruise weight to gross weight. 

^ 6w a 

W-TILDE = W = 

pu2 x]{l - yf)5/2 

6w 

Total Area = — 


I SLOPE = Code for signature type 

0 = Flat-top minimum overpressure 

1 = Minimum shock 

H, B, C, D, LAM, and yf are coefficients in equations for F-function and 

area distribution. (y^- = T; X = LAM) 

X = Axial position for F-function and area distribution near 

aircraft. 

X-BR =x-6z= Shifted axial position for pressure signature on the 

ground. 

DELTA P = Ap = Overpressure levels corresponding to x - Bz. 

P-j , P 2 = Overpressure levels in pressure signature. Pi and P 2 at 

same X-BR indicate a shock. 
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APPENDIX 


SAMPLE C;^E 


aHPUT ri=c 


:=b 0 0 0 0 , 0 1« UG=6 0 0 0 0 fi . n ;i L=3 0 0 - 0 I S= 0 ? I S: I G= 1 ? PERCEH= « 5 ? ‘!TL= . 1 


0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 GOOD 0 0 0 0 0 0 GOOD 0 0 0 0 0 0 GOOD 0 DO 000 

1 7 2 A 5 6 7 8 9 13 inM3 H 15 16 17 U I j -C '.i 27 7Z 74 22 76 27 7i 29 30 31 3? 33 34 35 36 37 38 3S 40 4' i7 43 44 45 46 47 48 43 50 51 5? 53 54 55 5S 57 58 55 60 61 62 63 64 53 63 67 64 t1 TO 7t 7? 73 74 75 76 77 73 79 80 

I n 1 11 11 1 1 n 1 11 1 1 11 11 1 1 1 i n ] n ] 1 1 n M ] 1 1 n ] 1 1 1 1 1 1 M 1 1 1 1 1 1 n 1 m 1 1 1 n i n ii n i n i 
22222222222 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 22222 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 
3 3333 33333 3 33333333 3 3333333333 3 3 3 33 3 33333 3333333 3333 333 3 33 3 3 

44444 444 4444444444444444444444444444444444444444444444444444444444444444444444 

555 555555555555555555555555555555555555555555555555555555555555 55 55 55555555 

6066666666 666666 6666666G 6 6666666666 666666666 G66G66B 668666666 8655 6 66G 

nil 111 inn nnnnnnn innnnnnnnnnnnnn nn innnin nn n 
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minimum-overpressure 


case an "IS” value of zero was used to produce both a 
solution and a minimum-shock solution. 
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(LAd-YF ) / (L-YF ) « • 820 95 9o It +00 

LOCATIOM OF FRONT SHOCK 16109 3 60 6^0 3 
LOCATION OF RtAR SHOCK . 5 A3 71 628 E +0 3 

TIME dtTWEEN SHOCKS .2 6 96 A908 c +00 

characteristic overpressure •21553217E+01 
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Figure 3.- The low-boom, high-drag paradox. 








Figure Special cases shown as function of weight factor 
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Figure 7.- Variation of shock and impulse with nose length yf. M = 2.7; h = 18 288 m (60 000 ft); 

W = 272 155 kg (600 000 lb); Z = 91.44 ra (300 ft). 
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Figure 8.- Estimated drag increments associated with configuration changes for sonic-boom minimization. 
M = 2.7; h = 18 288 m (60 000 ft); W = 272 155 kg (600 000 lb); 1 = 91.44 m (300 ft). 
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Figure 9.- Overpressure as a function of aircraft length. M = 2.7; W = 272 155 kg (600 000 lb); 

h = 1 8 288 m (60 000 ft) . 





1. R«porI No, 

NASA TP-1348 

2. Government Accession No. 

3. Recipient's Catalog No. 

4. Title and Subtitle 

SONIC-BOOM MINIMIZATION WITH NOSE-BLUNTNESS RELAXATION 

5, Report Date 

January 1979 

6. Performing Organization Code 

7. Author(s) 

Christine M. Darden 

8. Performing Organization Report No. 
L-12464 

10. Work Unit No. 

505-09-23-11 

9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, VA 23665 

11, Ck>ntract or Grant No. 

13. Type of Report and Period (Covered 
Technical Paper 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 

14. Sponsoring Agency Code j 

— 1 


15. Supplementary Notes 


16. Abstract 


A procedure which provides sonic-boom-minimizing equivalent area distributions 
for supersonic cruise conditions is described. This work extends previous 
analyses to permit relaxation of the extreme bluntness required by conventional 
low-boom shapes and includes propagation in a real atmosphere. The procedure 
provides area distributions which minimize either shock strength or overpressure. 


17. Key Words {Suggested by Author(s)) 

Sonic-boom minimization 
Impulse 

Pressure signature 
Equivalent area 


18. Distribution Statement 

Unclassified - Unlimited 


Subject Category 02 


19. Security Classif, (of this report) 
Unclassified 


20. Security Classif. (of this page) 
Unclassified 


21. No. of Pages 

50 


22. Price* 

$4.50 


^ For sale by the National Technical Information Service. Springfield. Virginia 2216! 


NASA-Langley , 1979 








